Lattice QCD determination of patterns of excited baryon states 
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Energies for excited isospin I ~ \ and I — § states that include the nucleon and A families of 
baryons are computed using quenched, anisotropic lattices. Baryon interpolating field operators that 
are used include nonlocal operators that provide G2 irreducible representations of the octahedral 
group. The decomposition of spin | or higher spin states is realized for the first time in a lattice 
QCD calculation. We observe patterns of degenerate energies in the irreducible representations of 
the octahedral group that correspond to the subduction of the continuum spin | or higher. The 
overall pattern of low-lying excited states corresponds well to the pattern of physical states subduced 
to the irreducible representations of the octahedral group. 
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I. INTRODUCTION 

The theoretical determination of the spectrum of 
baryon resonances from the fundamental quark and gluon 
degrees of freedom is an important goal for lattice QCD. 
To date there have been many lattice studies of ground 
state energies for different baryons [1-3] but only a few re- 
sults for excited state energies have been reported. [4-11] 
No clear determination of states with spin | or higher has 
been published because nonlocal operators have not been 
used. In this work we find degenerate energies that oc- 
cur in irreducible representations of the octahedral group 
corresponding to the subduction of the continuum spin 
I or higher. The pattern of lattice QCD baryon states 
that we observe is compared with the physical spectrum. 

Lattice correlation functions correspond to definite ir- 
reducible representations (irreps) of the octahedral group 
when the source and sink operators transform accord- 
ingly. There are six double-valued irreps of the oc- 
tahedral group: three for even-parity that are labeled 
with a g subscript (gerade) and three for odd-parity that 
are labeled with a u subscript (ungerade). They are: 
Gi g , H g , G2 g , G\ u , H u and G2 U - 

Continuum values of total angular momenta are real- 
ized in lattice simulations by patterns of degenerate en- 
ergies in the continuum limit that match the patterns in 
Table I for the subduction of spin J to the double- valued 
irreps of the octahedral group. For example, a state in 
one of the G2 irreps is a signal for the subduction of 
continuum spin | or higher. For spin |, there must be 
partner states in the H and G2 irreps that would be de- 
generate in the continuum limit. For spin |, there must 
be partner states in the G\, H and G2 irreps. 



TABLE I: The number of occurrences of double-valued irrep 
A of the octahedral group for half-integer values of continuum 
spin J. 
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This paper reports on work to determine the pattern 
of low- lying states in the I — \ and / = | channels. 
We carry out an analysis in quenched lattice QCD us- 
ing quasi-local and one-link-displaced operators formed 
from three quark fields. Smeared quark and gluon fields 
are used. [12, 13] Smearing reduces the couplings to 
short wavelength fluctuations of the theory and provides 
cleaner determinations of energies. [14] This is important 
when a large array of interpolating field operators is 
used in order to implement the variational method of 
Refs. [15, 16]. 

Both quasi-local and nonlocal operators transforming 
according to the G2 irrep are required. We have devel- 
oped sets of baryon operators that transform according 
to irreducible representations in two publications, one us- 
ing an analytical method based on appropriate Clebsch- 
Gordan coefficients for the octahedral group [17] and 
the other using a projection method that has been au- 
tomated and provides very large sets of operators [18]. 
Results obtained from both methods arc found to agree. 
In this work we use, for the positive parity channels, the 
three-quark operators defined in Tables VI, VII, and X 
of Ref. [17]. These compose a complete set of quasi-local 
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operators plus the simplest set of nonlocal operators hav- 
ing one quark displaced relative to the other two. 

The level ordering of baryon masses has been modelled 
using the spin-flavor SU(6) quark model. Negative par- 
ity multiplets are degenerate under SU(6) symmetry but 
the symmetry is broken by addition of a spin-spin contact 
interaction between quarks [19]. As a result, the lowest 
N(J P = i ) and iV(| ) states become less massive than 
the unperturbed levels and other states, such as iV(| ), 

(A^ ), and A(| ), become more massive. This pattern 
is observed in our lattice calculations. Degeneracies of 
SU(6) are further broken by addition of a spin-spin ten- 
sor interaction. In the positive parity excited states, Is- 
gur and Karl [20] introduced an anharmonic perturbation 
in order to explain the mass splittings, and a hypcrfinc 
interaction was introduced to break the degeneracy of the 
multiplets. They found that the second excited iV(i + ) 

state, the lowest 7V(| + ) state and the iV(| + ) states are 
nearly degenerate. Our lattice results using a relatively 
large pion mass equal to 490 MeV follow this pattern. 

This paper is organized as follows. In section II, we 
review the baryon operator construction and show the 
types of operators used in our simulations. In section III 
the use of charge conjugation together with time rever- 
sal is discussed. The combination provides a means to 
obtain a second correlation matrix for each gauge config- 
uration. In section IV we review the computational tech- 
niques with special emphasis on the variational method 
that is essential for obtaining excited states. Also we 
discuss the anisotropic action that is employed and the 
tuning of parameters required to obtain the correct speed 
of light. Section V presents and discusses the results for 
energies and eigenvectors in each irrep and their stabil- 
ity with respect to the number of operators used in the 
variational method. Section VI discusses the pattern of 
lowest energy levels in the lattice results and shows that 
it is consistent with the physical spectrum. Two notable 
findings are: 1.) the first positive-parity excited state 
in the I = \,G\ g channel is found to have energy sig- 
nificantly higher than the experimentally-known mass of 
the Roper resonance, and 2.) spin assignments of cither 
| or | are possible for some lattice levels that preferably 
would be assigned spin | in order to match the pattern 
of the physical levels. Section VII gives a summary of 
the results. 



II. IMPROVED BARYON OPERATORS 

A. Quasi-local operators 

We employ operators that transform as irreducible rep- 
resentations of the octahedral group. They are formed as 
linear combinations of elemental operators that have def- 
inite isospin and strangeness as follows, 

B fc AA) ( x .*) = C^fiB^yi^t), (1) 



TABLE II: Baryons and the corresponding quasi-local three- 
quark elemental operators. Columns 1 to 4 show the symbol, 
isospin, strangeness, and the form of elemental quasi-local 
source operators with maximum I z , respectively. In the next 
column, the MA label (mixed-antisymmetric) denotes com- 
binations of three Dirac indices that are antisymmetric un- 
der permutation of the first two labels and orthogonal to the 
totally antisymmetric combination. Similarly, the MS label 
(mixed-symmetric) denotes combinations of indices that are 
symmetric with respect to the first two labels and orthogo- 
nal to the totally symmetric combination. Columns 6-8 show 
numbers of embeddings of operators with irreps Gi g , Gig and 
H g , respectively. The last row shows the total number of 
even-parity operators, including all rows. Odd-parity opera- 
tors are obtained from the even-parity ones as explained in 
the text. 
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where a, /?, 7 are Dirac indices and the Clebsch-Gordan 
coefficients for the octahedral group provide the appro- 
priate coefficients c^^\. [17] The subscript '0' denotes 
quasi-local operators in which all quark fields have the 
same coordinates, (x, t). Octahedral group irreps are 
denoted by A = {G lg , G 2g , H g , G lu , G 2u , H u } and the 
corresponding irrep dimensions are d\ = {2, 2, 4, 2, 2, 4}, 
respectively. Each baryon operator carries a row label 
that distinguishes between the c?a members of irrep A, 
i.e., A = 1, 2, • • • , d,A- When group representations con- 
tain m\ occurrences of irrep A, the label k is used to 
denote the k-th embedding. 

Quasi-local elemental operators £? Q /3 7 in Eq. (1) are 
color singlet combinations of products of three quark 
fields, each of which is smeared in the same gauge- 
covariant fashion about point x. They are listed in Ta- 
ble II for the different baryons. Although not shown in 
the table, each quark field has a color index i, j or k, and 
a factor ajk is included in order to construct a color- 
singlet combination. The Dirac index symmetries that 
give non-vanishing operators and the number of avail- 
able quasi-local operators distributed over gerade irreps 
are also shown in the table. 

Note that the "barred" form of the operator is given 
and it is formed from three "barred" quark fields. The 
"unbarred" form of baryon operator is obtained by a sim- 
ilar sum using complex conjugates of the same coefficients 
and three quark fields, 

4 AA W)=ci A / 7 ?iWx>i)- (2) 
Note that no G 2 irreps are obtained from quasi-local in- 
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terpolating fields. 

A quark field operator and a "barred" field operator 
arc transformed into one another by charge conjugation 
as follows, 



Cq a tf = -C, 



act' Ha' i 



Cq a tf = q a ,C\, c 



(3) 



where C = 7472- A three-quark operator that transforms 
according to irrep A of the octahedral group is related 
by charge conjugation to an operator that transforms ac- 
cording to irrep A c , which has opposite parity, i.e., 



_ _ (AXk) p, R 

cBi A ^c^ci A ;;>cB a3 ^ 



(A c A c fe) 



„(AAfc) 
"a/37 



D °-'P'r Kj a'B''i'a 



JA c X c k)» 



(4) 



where C Q ' / g/ 7 ' Q( g 7 — Ca'aCpipCyj. If A is G\ g ,Gi g 
or then A c is G\ U ,G2 U or respectively, and 
A c = d,A + 1 — A. This one-to-one correspondence be- 
tween an operator of one parity and an operator with 
the opposite parity holds for any quasi-local or displaced 
operator. Choosing the gerade operators according to 
Eq. (1) and the ungerade operators according to the right 
side of Eq. (4) provides a convenient labeling of operators 
that facilitates the use of charge-conjugation relations to 
realize improved statistics, as explained in Sec. III. 



B. One-link-displaced operators 

A simple extension of quasi-local operators involves 
displacement of the third quark field along a spatial di- 
rection relative to the first two quark fields, with a gauge 
link included in order to maintain gauge covariance. An 
abbreviated notation for displacements uses the operator 
di whose action on a three-quark quasi-local operator is 
defined as follows, 



TABLE III: The number of occurrences of single-valued irrep 
A of the octahedral group for integer values of continuum 
angular momentum L. 
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follows, 

{A{B D+B dJb D B E B E 2 ~B) = 
I -±= (d x B + d y B + d z B + d- x B + d-yB + d-zB) \ 

§ [(d X B_- d_ x Bj + l {dyB- d-yB)] 

-i[(d x B - d- x B)- i{d y B-d_yB)} 
~^{d z B-d. z B) 
^[2{d z B + d- z B) - (d x B + d- x B) - {d y B + d- y B)] 



\ 



;[{d X B + d- X B) - {dyB + d-yB)] 



I 

(6) 



where displacement operators denoted as Ai,D± t o, and 
-^0,2 create combinations of operators that transform ac- 
cording to the A\ , T\ , and E single- valued irreps of the 
octahedral group, respectively. As may be seen in Ta- 
ble III, the lowest angular momentum in A\ is L = 0, the 
lowest in T\ is L = 1 and the lowest in E is L = 2. The 
combinations of displacements are chosen so that they 
transform as the lattice discretizations of the spherical 
harmonics Y Lm , i.e., A x ~ Too, £>+,o- ~ *n, *io, 
and -Eo,2 ~ Y20, (^22 + ^2-2)- Using this convention, the 
A\, Ti, and E one- link-displaced operators are defined 
as 



A r A 
A E A' 
A r A' 



P R (A ' V) 



_ (AXk) 



a{3~y 



(7) 



d e B ( t X \^t)=e t3k f a 



x t)qf{*, t% "(x + £, t)Uj K K (x, t), (5) 

where I can take any of six spatial directions, ±x, ±y, ±z. 

Displacement operators transform amongst themselves 
under lattice rotations. Linear combinations of displace- 
ments provide suitable bases for the construction of ir- 
reps. We define six linearly independent combinations 
of displacements and the operators that create them as 



(IS) (AXk) 
be apy 



ffe'fc/ 
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respectively, where C , I are Clebsch-Gordan 

VA A' X"J 

coefficients appropriate for forming overall irrep A from 
direct products of irreps A' and A". The coefficients 
follow the conventions of Ref. [17]. The second line of 
Eq. (7) defines new coefficients for linear combinations 
of displaced three-quark operators that realize the A\, 
Ti, or E-type operators of Eq. (6). Table IV shows 
the numbers of positive-parity, onc-link-displaced oper- 
ators for different baryons. For the one-link-displaced 
operators, there are inequivalent operators for each of 
the two ways that I — \ can be formed from three 
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TABLE IV: Numbers of embeddings of even-parity irreps, 
Gig, G2 g , and H g , that are obtained with one- link-displaced 
operators for different baryons. Columns show the number of 
operators for each overall irrep that can be made using the 
three irreps of one-link displacements: Ai, Ti, and E. Op- 
erators corresponding to a total derivative acting on a three- 
quark operator are omitted because they vanish when pro- 
jected to zero momentum. 
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light quarks, namely, the mixed- antisymmetric (MA) and 
mixed-symmetric (MS) combinations. Operators corre- 
sponding to total derivatives are excluded because they 
vanish when projected to zero total momentum. 

We construct matrices of correlation functions 



(AXk)* JA'A'fc'). 



Cj£?'(t)5AA'<SAA' 

(0,0)10)7*^7^,7^, 



(8) 



where displacement operator de applied to baryon 
fields denotes a one-link-displaced baryon field for £ = 
±x, ±y, ±z and a quasi-local baryon for I = 0. Lattice 
operators belonging to different irreps or different rows 
are orthogonal because of the octahedral symmetry of 
the lattice. Different embeddings of a given irrep and 
row provide sets of operators that we used in calculat- 
ing the matrix of correlation functions. The 74 matrices 
are included in Eq. (8) in order to produce a Hermitian 
matrix of correlation functions. For the operators used 
in this work, the 74 matrices always reduce to a factor 
Vffl — ±1 when they act on the three "barred" quark 
fields in B (AX) 



k . This factor is the p-parity defined in 
Ref. [17]. We write a matrix of correlation function for 
irrep A and row A as 



>(AA) 



(x,t)B^ A) (O,0)|0). 



(9) 



III. CHARGE CONJUGATION 

By inserting complete sets of particle states \n) and 
antiparticle states |n) in Eq. (9), and using translational 
invariance to extract the dependence on x and t, corre- 
lation functions are expanded as 

CjfrHt) = P fe ^P,o[^^(t)(0|4 AA V)(n|^ A) |0)e- B " t 

n 

-^^(-0(0|4' A V)<n|4 AA) |0>e^]. (10) 

n 

Using the charge conjugation relations \n) = C|n)e l< ^, 
and the invariance of the vacuum state under charge con- 
jugation, the antiparticle contributions in the t < part 
of Eq. (10) are rewritten as 

V ( k K) {0\CB { v X) tfC\n){n\tfCB ( k K A) C t |0> 



where 



The relation between correlation 



functions with different parities is then [18], 

-,(A C A C )*/ 



(11) 



Applying the temporal lattice boundary conditions 

clt A) M) = VtCii X) {T - t), with T) t = +1 for periodic 
boundary conditions and rj t — — 1 for antiperiodic bound- 
ary conditions, the matrix of correlation functions can be 
written as 



a 



(AA), 



^ ; (*)=^p,oEK*)H' ) <0|B 

71 

-ry t 0(T-O^ c ^Ol4 VcAc V}>l^ cAc) |Ore~^ (T ~ t) l.(12) 



in the interval < t < T. The forward propagating signal 
of a correlation function is equal to the backward propa- 
gating signal of the parity-reversed, complex-conjugated 
correlation function within the factor — t] t , i.e., 



t). 



(13) 



This symmetry is used to improve statistics. For each 
gauge configuration we compute matrices of correlation 

functions using for positive-parity operators the B^- 
type operators and for the corresponding negative-parity 
operators the B^"^ -type operators. Starting with a ma- 
trix of correlation functions for a given parity, a second 
matrix of correlation functions is obtained for the oppo- 
site parity. After time-reversal, complex conjugation and 
multiplication by the factor —Tit, the second matrix of 
correlation functions is averaged with the first one. A 

typical state created by the B k AX ^ operator has an ef- 
fective mass plateau occurring in first half of the time 
extent, t < T/2. Using the B k A " Xc ^ operator, a simi- 
lar plateau occurs at t > T/2. These two plateaus are 
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largely statistically independent samples because gauge 
configurations on time slices significantly before T/2 are 
largely uncorrelated with ones on time slices significantly 
after T/2. 



IV. COMPUTATIONAL TECHNIQUES 
A. The variational method 

Analysis of excited state energies is based upon the ma- 
trices of correlation functions defined in Eq. (9). Statis- 
tics are improved by averaging over the rows of irreps, all 
of which are equivalent because of the cubic symmetry. 
The matrix averaged over rows is denoted as cj£.} (t), i.e., 
the row label is omitted. 

Because the operators used to form a matrix of corre- 
lation functions are not normalized with respect to one 
another, we find that they can produce diagonal elements 
of a correlation matrix that differ by two orders of mag- 
nitude. It is convenient to adopt a normalization scheme 
at time to such that each operator produces a diagonal 
matrix element equal to 1. Thus, we define normalization 
factors, 



N k (t ) 



(14) 



and renormalize the operators by attaching a factor N k 
to operator B k . The new matrix of correlation functions 



is 



C$(t)=N k C$Wk> 



(15) 



and its diagonal elements obey C^\to) = 1. Because of 
this normalization convention the components of eigen- 
vectors indicate the relative importance of the contribu- 
tions of various operators to an eigenstate. 

In order to extract an energy spectrum from the matrix 
of correlation functions, we numerically solve the follow- 
ing generalized eigenvalue equation, 

£^(t)4?>(t,t ) - a^\tMYAihto)v£\t,t Q ), 
h 



/ j ^kk< 

k> 



(16) 

where superscript n labels the eigenstates, such as the 
ground state, the first excited state, and so forth. The 
reference time to m Eq. (16) is taken near the source time 
t = in order to have significant contributions from ex- 
cited states and to ensure stability of the Cholesky pro- 
cedure below. The principal eigenvalues a^ n '(t,to) are 
related to the energy E n by [16] 

aW (t, t Q ) ~ e - E ~ (t - to) (l + 0(e-l 5£ l*)) , (17) 

where SE is the difference between E n and the next clos- 
est energy. 



Energies E n are calculated from the generalized eigen- 
values according to 



E„ = - In 



«<")(* + Mo) 



(18) 



We also perform fits of cr") (i, to) to an exponential func- 
tion a( n \t, to) = e~ En ^~ to ^ over a range of time slices in 
order to better determine the error of the energy E n . 

Numerical solutions of the generalized eigenvalue 
equation are obtained by performing a Cholesky- 
decomposition of the matrix of correlation functions at 
the reference time, to, 



C { k K k ){to) = A T kk ,,A k ,, k ,. 
For the transfer matrix, 

T(t,t ) = (A T )- 1 C (A) (t)A- 1 . 
the eigenvalue problem is 

T kk ,(t,to)v£r\t,to) = « w (Mo)n (n) (Mo), 
v£ n \t,t ) = A kk ,v$ ) (t,t ). 



(19) 

(20) 

(21) 
(22) 



Left and right eigenvectors of the transfer matrix 
T kk '(t, to) are the same. 

Eigenvectors in the generalized eigenvalue equation, 
Eq. (16), are orthogonal with respect to Cj^}(to), i-e., 

V k n)T (tMCli)(to)v<£'\t,to) = Snn'. (23) 

It follows from Eqs. (16) and (23) that the correlation 
matrix at time t is diagonalizcd by these vectors, i.e., 



S n ) T t 



(t,to)Citht)v { « ) (t,to) = aW(t,t )6< 



(n). 



(24) 



It is useful to define improved operators by first adopting 
a suitable normalization of the eigenvector components. 
Let 



v k n \t,to) = Z n {t)N k v k n, {t,t ), 



(«), 



such that 



5>^(Mo)| 2 = l. 



(25) 



(26) 



The factor Z n allows a comparison of the overall normal- 
ization of eigenvectors calculated on different volumes. 



B. Lattice setup and tuning 

A technical difficulty in extracting excited state ener- 
gies from lattice QCD simulations is that the signal-to- 
noise ratio degrades for increasing time. A plateau in the 
energy may not be realized for excitations. We use sev- 
eral techniques to overcome this problem, namely, the op- 
timization of operators so as to achieve early plateaus of 
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the energies, the imposition of symmetries and the use of 
more gauge configurations to increase the statistics, and 
the use of anisotropic lattices. Anisotropic lattices are de- 
signed with temporal lattice spacing a t smaller than spa- 
tial lattice spacing a s . Using a finer lattice spacing along 
the time direction provides more time-slices for analysis 
of correlation functions at small time separations from 
the source. Of course anisotropic lattices require tuning 
of coefficients in the action. 

In this work, anisotropic lattices with two different 
volumes are used: 239 gauge field configurations on a 
16 3 x 64 lattice and 167 configurations on a 24 3 x 64 lat- 
tice. For both lattices, the renormalized ratio of spatial 
lattice spacing to temporal spacing is £ = a s /a t — 3.0 and 
the temporal lattice spacing corresponds to a^ 1 = 6.0 
GeV [8]. Gauge-field configurations are generated using 
the anisotropic, unimproved Wilson gauge action in the 
quenched approximation with (3 — 6.1 



1. Anisotropic Wilson gauge action 



by 



The anisotropic Wilson gauge action [21-24] is given 



S G=£- c [l E RcTr(l-P ss ,(x)) 



x,s>s' 



•£o]TRcTr(l-.P ;st (x)) 



(27) 



The renormalized anisotropy £ = 3 is held fixed in our 
calculations and the bare anisotropy, £□> is varied in order 
to obtain the target value of £ for the desired value of (3 = 
6.1. Because the quenched approximation is used, tuning 
of the gauge action may be performed independently of 
the fermion action. Determination of the renormalized 
anisotropy is based on measurements of static potentials 
for a quark-antiquark pair. 



2. Anisotropic Wilson fermion action 
The anisotropic Wilson fermion action has the form, 
S| = a t a 3 s ^2~q(x) [m + v t (74V* - y A t) 

X 

+v s J2 (7.V. - ^A 8 )] (28) 

s 

where 

a , , U li (x)q(x + p,) + UUx-p,)q(x-p,)-2q(x) 

A* 

U„(x)q(x + (i) - W(x - (i)q(x - fi) 
V^(x) = ^ . (29) 

Note that when v t = v s , the fermion action in Eq. (28) 
becomes the original Wilson fermion action [25] . We hold 




FIG. 1: Pion dispersion relation based on tuned parameters 
mo and v a . The renormalized anisotropy is 3.0 and the lattice 
volume is 16 3 x 64 with [3 — 6.1. Lines show the the c = 
1 dispersion relation. The meaning of symbols is given in 
Table V. 



u t = 1 in the tuning because that preserves the projection 
property of (1 ± 7 4 )/2 for parity. Then v 8 and the bare 
quark mass mo arc tuned in order to obtain the desired 
pion mass and the correct speed of light. 

The tuning is based on the relativistic dispersion rela- 
tion for the pion, E 2 (p) = m 2 c 4 + c 2 p 2 , where E(p) 
is the energy of a pion with total momentum p and 
mc 2 = E(0). Expressing quantities in terms of the lattice 
spacing, the speed of light is 



c( P ) 



' (a t E(p)) 2 - (a t E(0)y 



(a s p) 2 



(30) 



where £ is the renormalized anisotropy. We measured the 
pion energy E(p) for p = and the lowest three nonzero 
values of momentum, namely, p = -^(1,0,0), -^(1, 1,0), 



and 



^7(1, 1, 1). We tuned v s for five different bare quark 
masses mo- The dispersion relation for tuned parameter 
v s on the 16 3 x 64 lattice is plotted in Fig. 1. The straight 
lines passing through the mean value of E 2 (0) show the 
desired continuum dispersion relation for each pion mass. 
In physical units, the volume is about (1.6 fm) 3 x (2.1 fm). 
For the anisotropic, unimproved Wilson fermion action, 



TABLE V: Tuning parameters for the fermion action on the 
16 3 x 64 lattice using renormalized anisotropy £ = 3 and v t = 1 



Symbol 


2 2 

a t m n 


m w (MeV) 


m 







0.0066(5) 


490 


-0.313 


0.898 


□ 


0.0104(5) 


610 


-0.318 


0.901 


A 


0.0126(5) 


675 


-0.305 


0.902 


O 


0.0224(5) 


900 


-0.290 


0.910 


< 


0.0335(5) 


1100 


-0.280 


0.915 



tuning of the fermion action produced the results given 
in Table V. A periodic boundary condition is employed 
for spatial directions and an anti-periodic boundary con- 
dition is employed for the temporal direction. Our spec- 



7 



trum calculations are based on the parameter set that 
yields a pion mass of 490 MeV. 

Gauge links are smeared at source and sink according 
to the APE-smearing method [12], 

U i n \ x ) + I E Ui n) (x)uj n \x + k)ul {n \x+j), (31) 

where projection of gauge link matrices to SU(3) is per- 
formed after each iteration. We used the APE-smearing 
parameters (a,n) = (2.5,3) throughout this work. 

Quark fields are smeared for both quasi-local and one- 
link-displaced operators using the Gaussian smearing 
method [13], 

qil (x,t) ^Y>G^\x,a/) qil ^,t), (32) 

x> 

where is an operator that acts recursively, 

& N \x,x') = Y,{^.,y + ^l, y /^N)G {N - 1 \y,x'), 
y 

GW(x,x') = 6 x>x ,, (33) 

and x , is a three-dimensional, gauge-covariant Lapla- 
cian operator. Wc used (a, N) = (3.0, 20) in this calcu- 
lation. 

V. DISCUSSION OF RESULTS FOR EACH 
SYMMETRY CHANNEL 

We have extracted energies for I = i and I = | 
channels by diagonalizing matrices of correlation func- 
tions formed from three-quark operators that share the 
same octahedral symmetry. A large number of operators 
is available so that reduction to a set of the most im- 
portant operators is carried out in the initial stage. The 
available operators are divided into several subsets con- 
taining roughly 3 to 8 operators each and we diagonalize 
the matrices computed from each subset of operators. In- 
spection of the energies and corresponding eigenvectors 
reveals which operators provide stronger signals for the 
low- lying states of interest. In the second stage we omit 
all but these good operators and form new matrices of 
correlation functions, which are then diagonalized. This 
two-stage procedure yields solid results for energies of 
low-lying states and is less susceptible to noise than ap- 
plying the variational method directly to the matrices of 
largest dimension. 

In each symmetry channel, we have studied how the 
energies depend upon the number of good operators that 
are used in applying the variational method. Starting 
from a given set of good operators, matrices of corre- 
lation functions of reduced dimension are obtained by 
omitting a less important operator from the selected set. 



Continuing this process, we have diagonalized matrices of 
various dimensions and obtained energies. Low- lying en- 
ergies in each symmetry channel are plotted as a function 
of matrix dimension in Fig. 2. The I = |, G 2g / U channel 
is not shown because results in that channel are based 
on a single operator. The lowest energy states are more 
or less stable as long as a few important operators arc 
retained. However, it is expected that use of more than 
the minimum number of operators improves the lowest- 
lying states because it reduces contamination from higher 
lying states. This expectation appears to hold true for 
the lowest and first excited states but is less evident for 
higher excited states with the limited set of operators 
used in this work. 

Diagonalization of a matrix of correlation functions 
provides the principal eigenvalues a^ n \t, t ). These have 
been fit with a single exponential form over a selected 
range of times, t\ < t < t 2 in order to extract energies, 
i.e., a( n \t,t ) = Ce~ En ^~ to \ The resulting energies are 
given in Tables VI, VII , VIII, and IX, which show the 
energies obtained for positive-parity states with I = ^, 
negative-parity states with I = |, positive-parity states 
with / = | , and negative-parity states with I = |, re- 
spectively. Operators are selected so as to provide the 
best quality of energies for the larger lattice volume and 
the same sets of operators are then used for the smaller 
lattice volume. The Z-factors given in Eq. (25) are also 
shown along with the energies. Each table shows results 
for two lattices, 16 3 x 64 and 24 3 x 64, and six symme- 
try channels, G\ g / U , G2 g / U and H g / U . In some channels, 
there are states in the 16 3 x 64 analysis that are missing 
from the 24 3 x 64 analysis and vice versa. 

We have examined the Z-factors defined in Eqs. (25) 
and (26) as a possible means to detect scattering states 
following the methods of Ref. [26] . This has not produced 
clearly interprctable results. Therefore we discuss the re- 
sults for each channel in the following subsections assum- 
ing that there are no scattering states in the extracted 
spectra for the pion mass used in these calculations. 

A. The Gi channel 

In the I = \ , Gi g channel, we selected 10 operators on 
the 24 3 x 64 lattice that give the best quality of ener- 
gies using the variational method, and the same 10 op- 
erators are used to analyze the 16 3 x 64 lattice data. 
The ground state of this channel, i.e., the nucleon state, 
has a very stable set of eigenvector components with re- 
spect to time as shown in Fig. 3. As many studies have 
shown, the ground state is dominated by a local operator 
of the form, u 7 '(C^dP+u, where P± are parity projec- 
tion matrices for Dirac spinors. However, we find non- 
ncgligiblc contributions from nonlocal T\ operators repre- 
sented by diamonds, circles and inverted triangles in the 
figure. One significant operator has a mixed-symmetric 
isospin symmetry with totally symmetric Dirac indices. 
Although this isospin symmetry is unusual, nonlocal op- 
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FIG. 2: Energies in lattice units vs. dimension of the matrices of correlation functions. Data for the 24 3 x 64 lattice are used. 
Upper panels show energies for positive parity states and lower panels show energies for negative parity states. The / = | , G2 g / U 
states are excluded because correlation functions are based on a single operator. 
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FIG. 3: Eigenvector components v ® (t,0) of the ground state for (a) the I = |, Gi 9 channel using a 10 x 10 matrix and (b) the 
I = i, H u channel using an 8 x 8 matrix, both on the 24 3 x 64 lattice. The vertical axis is the value of the coefficient of a basis 
operator and the horizontal axis is time. Each basis operator corresponds to a symbol in the plot; the four most important 
operators are listed below the plot together with their symbols. Ns t s z is om " notation for a smeared three-quark operator with 
spin S, spin projection S z , irrep A, and embedding k. The symmetry (MA, MS or S) of Dirac indices for each three-quark 

operator is written next to the operator. Explicit forms of the operators Ns t s z are provided in Ref. [17]. 



erators formed from three quarks allow such combina- 
tions. In fact, this operator plays a significant role in the 
first excited state in combination with a local operator 
of the second embedding. The first and second excited 
states have relatively clean signals and are dominated 
by local and Ti-displaced operators, respectively. Sim- 



ilar energies are calculated for states on the two lattice 
volumes except for the second excited state. Dominance 
of a nonlocal operator in the second-excited-state eigen- 
vector might explain this dependence on lattice volume 
because we find that the nonlocal operators generally 
exhibit more sensitivity to the spatial volume. Eigen- 
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TABLE VI: The I = \ N* , positive-parity energies in the 
units of aT" 1 obtained on a 24 3 x 64 lattice (left) and a 16 3 x 64 
lattice (right). The irrep is shown in the first row together 
with the dimension of the matrix of correlation functions used. 
In the first column, energy values obtained from fitting the 
generalized eigenvalues are shown in increasing order, namely, 
the second row contains energies of the lowest-lying states for 
each channel and the third row contains the first excited state 
energies, etc. The Z-factor of Eq. (25) is given in the second 
column. A time range in which the energy and the Z-factor 
arc fitted on the 24 x 64 lattice is listed in the third column. 

Gi g , 10 x 10, 24 3 x 64 G lg , 10 x 10, 16 3 x 64 

energy Z n time energy Z n 

0.193(3) 5.75(15) xl0~ 5 19-24 0.194(4) 8.45(15) xl0~ 5 
0.398(6) 4.79(12) xl0~ 5 8-12 0.403(6) 9.53(26) xl0~ 5 
0.399(7) 5.64(33) xl0~ 5 10-15 0.409(8) 8.84(46) xl0~ 5 

G 2g , 4 x 4, 24 3 x 64 G 2g , 4x4, 16 3 x 64 

energy Z n time energy Z n 

0.377(13) 1.47(57) xl0~ 5 13-18 0.414(13) 2.11(35) x 10" 6 
0.398(10) 1.34(25) xl0~ 5 10-14 0.444(7) 7.66(100) x 10" 5 
0.574(2) 3.73(30) xl0~ 5 11-14 0.564(12) 5.90(199) x 10 -5 

H g , 12 x 12, 24 3 x 64 H g , 12 x 12, 16 3 x 64 

energy Z„ time energy Z n 

0.365(9) 5.40(15) xl0~ 5 13-18 0.375(9) 9.25(33) xl0~ 5 
0.395(7) 5.14(14) xl0~ 5 10-15 0.407(5) 8.38(21) xl0~ 5 
0.400(7) 4.39(7) xl0~ 5 10-16 0.407(8) 7.75(10) xl0~ 5 
0.403(8) 4.27(10) xl0~ 5 10-14 0.422(6) 8.07(17) xl0~ 5 

0.427(8) 6.36(8) xl0~ 5 
0.443(6) 8.75(7) xl0~ 5 



vectors are more or less the same in each state on the 
two lattices as far as the coefficients of the two or three 
most important operators are concerned. Fluctuations of 
the eigenvector components generally increase with time 
for excited states, but they are fairly stable in the time 
ranges where energies are extracted. 

In the G\ u channel we obtain two states which could 
correspond to AT(±~,1535) and AT(|~, 1650). Each is 
dominated by two local operators but in different linear 
combinations. Contributions from nonlocal operators for 
the first excited state are around \v^\t, 0)| 0.3, which 
is less than for positive parity states. The energy of the 
lowest state on the two lattce volumes is essentially the 
same, suggesting that a 16 3 lattice is large enough to 
contain this state with a 490 MeV pion mass. 

In the I — |,Gi g channel, the first excited (the low- 
est) state in 16 3 volume corresponds to the lowest (the 
first excited) state in 24 3 volume. This correspondence is 
based on the similarity of eigenvector components. The 
first excited state of I = has lattice energy of 

0.392(7) on the larger volume and 0.373(8) on the smaller 
volume, with error bars that do not overlap. The appar- 
ent crossover of energies of the same state at different 
lattice volumes is observed only in this particular state. 



TABLE VII: The I = \ N* , negative-parity energies in the 
units of aT" 1 obtained on a 24 3 x 64 lattice (left) and a 16 3 x 64 
lattice (right). A similar description applies as in the caption 
of Table VI. 



Gin, 7 x 7, 24 3 x 64 




Gm, 7x7, 16 3 x 64 


energy Z„ 


time 


energy Z„ 


0.290(5) 5.68(53) xlO™ 5 


12-17 


0.290(7) 9.19(40) xl0~ 5 


0.323(12) 5.73(44) xl0~ 5 


11-17 


0.340(6) 9.35(73) xl0~ 5 


G 2u , 3 x 3, 24 3 x 64 




G 2u , 3x3, 16 3 x 64 


energy Z n 


time 


energy Z n 


0.327(10) 1.93(82) xl0~ 5 


15-20 


0.325(9) 2.67(121) xl0~ 5 


0.473(11) 7.24(22) xl0~ 6 


7-12 




H u , 8 x 8, 24 3 x 64 




H u , 8x8, 16 3 x 64 


energy Z„ 


time 


energy Z„ 


0.284(16) 5.00(38) xl0~ 5 


18-23 


0.284(9) 8.07(26) xl0~ 5 


0.308(7) 5.74(8) xl0~ 5 


15-22 


0.310(8) 1.04(1) xl0~ 4 


0.330(9) 7.93(5) xl0~ 5 


16-22 


0.328(8) 1.44(l)xl0" 6 



TABLE VIII: The I = § A*, positive-parity energies in the 
units of aT" obtained on the 24 x 64 lattice (left) and the 
16 3 x 64 lattice (right). A similar description applies as in the 
caption of Table VI. 



Gi g , 5 x 5, 24 3 x 64 


Gi 9 , 5x5, 16 3 x 64 


energy Z n 


time 


energy Z n 


0.378(17) 5.05(14) xl0~ 5 


14-18 


0.373(8) 1.67(5) xl0~ 4 


0.392(7) 7.41(26) xlO -5 


12-17 


0.401(13) 9.33(16) xl0~ 5 


G 2g , 1 x 1, 24 3 x 64 


G 2g , 1 x 1, 16 3 x 64 


energy Z n 


time 


energy Z n 


0.416(7) 8.92(28) xl0~ 5 


9-14 


0.429(6) 1.49(4) xl0~ 4 


H g , 4 x 4, 24 3 x 64 




H g , 4x4, 16 3 x 64 


energy Z„ 


time 


energy Z„ 


0.231(3) 1.22(3) xl0~ 4 


19-26 


0.232(4) 1.79(4) xl0~ 4 


0.326(7) 5.89(5) xl0~ 5 


13-20 


0.341(13) l.ll(l)xl0" 4 


0.370(25) 6.54(4) xl0~ 5 


19-23 


0.399(10) 1.20(6) xl0~ 4 


0.400(8) 7.15(5) xl0~ 5 


14-19 


0.402(6) 1.29(2) xl0~ 4 



The eigenvectors for these states show that they are dom- 
inated by T\ displaced operators with MS isospin. 

In the negative parity channel, it is more difficult to 
determine the excited state energies due to the large 
error bars, while the lowest state, corresponding to 
A(| ,1620), is fairly stable. Again the T\ displaced 
operators dominate the eigenvectors of the two lowest 
energy states. 



B. The H channel 

We find a richer spectrum in / = i , H g / U than in other 
channels as shown in Tables VI, VII, and VIII. Plateaus 
are achieved for t w 20 in the first three states. Approxi- 
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TABLE IX: The I = § A*, negative-parity energies in the 
units of a^ 1 obtained on the 24 3 x 64 lattice (left) and the 
16 3 x 64 lattice (right). A similar description applies as in the 
caption of Table VI. 



Gm, 6 x 6, 24 3 x 64 
energy Z n time 

0.309(9) 3.13(99) xl0~ 6 15-20 
G 2u , 1 x 1, 24 3 x 64 
energy Z n time 

0.506(7) 7.31(16) xl0~ 5 6-12 

H u , 7 x 7, 24 3 x 64 
energy Z n time 

0.320(5) 8.38(30) xl0~ 5 12-17 
0.485(7) 5.56(53) x!0~ 5 7-12 



Giu, 6x6, 16 3 x 64 
energy Z„ 
0.312(7) 5.38(117) xl0~ 5 
G 2u , 1 x 1, 16 3 x 64 
energy Z n 
0.539(5) 1.43(2) xl0~ 4 

H u , 7x7, 16 3 x 64 
energy Z n 
0.317(6) 1.16(4)xl0~ 4 
0.485(8) 9.14(28) xl0~ 5 



mate plateaus could be selected for the next three states. 
However, the use of more diverse operators together with 
better statistics is required in order to extract reliably 
more than the first three states. 

We obtained four energies for the I = \, H g channel 
using the 24 3 lattice and six energies using the 16 3 lat- 
tice, with errors relatively small compared with those in 
other channels. A quasi-local operator yields the dom- 
inant coupling to the lowest H g state and our opera- 
tor is essentially the same as the Rarita-Schwinger pro- 
jected operator that has been used in other works. The 
Rarita-Schwinger projected lattice interpolating field be- 
longs to the H irrep and is therefore consistent with 
the subduction of continuum spins |, |, |, • • • . Several 
works [4, 6, 7, 9] have suggested that this lowest H g 

state likely corresponds to a continuum spin | state. 
However, we find an indication that it corresponds to a 
higher spin, either | + as in the iV(| + , 1680), or | + . The 
possible spins will be discussed in the next section when 
the patterns of states over all the irreps are considered. 

For states above the ground state, combinations of 
nonlocal T\ operators dominate the couplings. Consis- 
tent with this, we observe a finite volume effect on some 
energies in the N* , H g excited spectrum. Even though 
the lowest energy state is dominated by the quasi-local 
operator, its energy has some dependence on the volume 
of the lattice. The lowest and the first excited states 
are clearly identifiable from eigenvector compositions in 
both volumes. Eigenvectors for the second and the third 
excited states are contaminated by noise on the 24 3 lat- 
tice with the available statistics, although the 16 3 lattice 
provides stable eigenvectors. 

The / = \ , H u states provide the least contaminated 
group of states in the sense that energies for several ex- 
cited states are determined easily and the lowest few 
states have eigenvectors that are nearly as stable as those 
of the nucleon G\ g ground state. The lowest-lying state 

should correspond to the physical 7V(| , 1520) state, the 



first excited state should correspond to the Af(| , 1675) 
state, and the second excited state should correspond 
to the A^d , 1700) state, supposing no scattering states 
are involved. Energy differences between these states are 
relatively small and our lattice results agree with this 
pattern. In contrast to the positive parity channel, finite 
volume effects in energies are well within statistical fluc- 
tuations. The nonlocal three-quark operators that we use 
yield nonvanishing MS isospin combinations although the 
simplest quark model does not have them. Our results 
show that contributions of MS isospin operators are cru- 
cial to the first and second H u excited states. The same 
operators are important in both volumes. Energy dif- 
ferences between parity-partner states are found to be 
smaller for the 1 =\,H irreps than for other channels. 

The A baryon appears as the lowest energy state in the 
I = | , H g channel. A local operator with nonrelativistic 
spin components dominates the ground state with very 
small errorbars. The first and second excited state sig- 
nals should correspond to A(| + , 1600) and A(| + , 1905), 
respectively. Signals for these states are also reasonably 
clean on the 16 3 lattice, however the results on the 24 3 
lattice are subject to larger fluctuations and the extracted 
energies are not reliable. 

In the / = \,H U spectrum we find that the ground 

state, corresponding to A(| ,1700), has stable signa- 
tures in both energy and eigenvector plots but the signal 
is short-lived. This is associated with the fact that the 
most important operator is the negative parity transform 
of the operator that dominates the positive parity ground 
state. The resulting long-lived plateau in the backward- 
in-time part of the correlation function limits the range of 
times where the excited state energies can be extracted. 
The energy differences between parity-partner states of 
the low-lying delta baryon spectrum in the H irrep are 
larger than in other channels. 



C. The G 2 channel 

The Gi baryon spectrum is the least explored because 
nonlocal operators have not been readily available. How- 
ever, G2 operators are very important for the assignment 
of higher spins because of their role in the patterns of 
subduction of continuum spins to the octahedral irreps. 

We have selected the four best operators for the / = 
\iGig spectrum by the method described earlier. One 
of these is the operator with .E-type spatial displace- 
ment, which transforms evenly under spatial inversion. 
We find that G2 energies are the most difficult to ex- 
tract in our calculations because the plateaus are rel- 
atively short-lived, signal-to-noise ratios are small, and 
the number of operators that we have used for the varia- 
tional method is limited. The limitation on the number 
of operators arises because we have restricted the oper- 
ators to ones that can be constructed from one-link dis- 
placements. Work in progress includes more varied types 
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of operators, which yield larger numbers of G 2 opera- 
tors. [11] 

The lowest state in the / = \-,G 2g channel has sig- 
nificant couplings to the Ti-displaced MA isospin oper- 
ator, the Ti-displaced MS isospin operator, and the E- 
displaced MS isospin operator. There is a non-negligible 
contribution to the lowest energy state of G 2 from the E- 
type of displaced operator that corresponds to L > 2 in 
the subduction of continuum angular momenta. The first 
excited state has significant couplings to a combination 
of two Ti operators that also contribute to the ground 
state but in a different linear combination. The second 
excited state couples almost purely to a single T\ opera- 
tor with MS isospin. Eigenvectors are similar on the two 
lattice volumes for the lowest three states, however their 
energies change significantly. The lowest energy is calcu- 
lated to be 0.362(17) in this channel, which agrees very 
well with the lowest I — \,H g energy, 0.365(9). Because 
these two energies are degenerate within errors, they are 
consistent with being the partner states for spin | . How- 
ever, that is not the only possibility as is discussed in the 
next section. 

For the I = \,G 2u channel, the Ti-displaced op- 
erator with MS isospin dominates the lowest state, 
while the _E-displaced operator plays a significant role 

with \v^\t : Q)\ ~ 0.4 in the first excited state. The 
lowest physical state corresponding to I = \,G 2u is 

7V(§ , 1675), which should also occur as the first excited 
state of the H u spectrum. The lowest energy state of 
H u should correspond to iV(| , 1520). We find that the 
first excited state of the H u channel has energy 0.308(7), 
whereas the lowest state G 2u state has energy 0.327(10) 
in the larger volume. In the smaller lattice volume, H u 
and G 2u energies are 0.310(8) and 0.325(9), respectively. 
The possible spin assignments are discussed further in 
the next section. 

Having several good operators in each symmetry chan- 
nel is important for the success of the variational method. 
Based on our present results using four G2 operators, we 
find that the energy of / = |, G 2u is lighter than that 
of I = i, G2 g , which is also consistent with experiment. 
An earlier analysis with one G 2 operator found the op- 
posite. [8] 

For the I — |, G 2 channels, our set of quasi-local and 
one-link-displaced operators does not contain a subset 
that provides better signals than are obtained from a cor- 
relation function based on a single operator. The problem 
is noise. Therefore we simply fit the energies based on in- 
dividual correlation functions and pick the best-behaved 
one, which is found to involve the Ti-displaced operator. 
Because no diagonalization is involved, the energy is rel- 
atively contaminated in early time-slices and the plateau 
is weak. The I — § , G 2 spectrum needs further study us- 
ing a larger set of operators, such as two-link-displaced 
operators or operators that have two quarks displaced 
from the third quark. 



VI. PATTERN OF LOWEST-LYING ENERGIES 

The physical spectrum of baryon excited states shows a 
number of degeneracies between states of different spins. 
In the lattice results, this means that particular care must 
be exercised in order to identify spins because acciden- 
tal degeneracies of two states can provide the same pat- 
terns in the octahedral irreps as a single higher spin state. 
Given the limited energy resolution of our calculations, 
this leads to alternative interpretations for the spins of 
some excited states in the lattice results. It is more ap- 
propriate to compare the pattern of lattice results to the 
pattern of experimental masses and spins subduced to the 
irreps of the octahedral group. Although there are sub- 
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FIG. 4: The energies obtained for each symmetry channel of 
/ = I baryons are shown based on the 24 3 x 64 lattice data. 
The scale on the left side shows energies in lattice units and 
the scale on the right side shows energies in GeV. The scale 
was set using the string tension. Errors are indicated by the 
vertical size of the box. 

stantial discretization errors with the quark action that 
is used, and they could contribute differently in the dif- 
ferent irreps, clear patterns in the degeneracies emerge. 
Focusing on the group of three positive-parity excited 
states near lattice energy atE — 0.36 in Fig. 4, two inter- 
pretations are possible. A.) The group consists of a spin 
h state and a spin | state that accidentally are degener- 
ate. In this case the G\ g state corresponds to spin | and 
the H g and G 2g partner states correspond to the subduc- 
tion of spin |. B.) The group consists of a single state 
with the degenerate G\ gi H g and G 2g partner states cor- 
responding to the subduction of spin |. Note that these 
are the only possible interpretations. It is not possible 
for the H state to be an isolated spin | state because 
that would require the G 2 state also to be an isolated 
state. There is no interpretation for an isolated G 2 state. 
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Higher spins than | would require more partner states 
than are found. 

In the physical spectrum of positive-parity nucleon res- 
onances, the lowest excited state, iV(i + , 1440), lies be- 
low all negative parity states. We do not find a signal for 
a positive-parity excitation that has lower energy than 
the negative-parity excitations at this quark mass. The 
next two excited nucleon states are essentially degener- 
ate, namely, iV(§ + , 1680) and A(± + , 1710), each with a 
width of about 100 MeV. Spin | states occur only at 
significantly higher energy (1990 MeV). These states are 
well separated from other states and the width of each 
is about 100 MeV. Primarily because of the absence of 
spin | in the low-lying spectrum, interpretation A.) of 
our lattice results is more consistent with the pattern of 
physical energies and spins. 

In the negative-parity spectrum, we also obtain essen- 
tially the same results for both lattice volumes. The 
three lowest states shown on the right half of Fig. 4 
are unambiguously identified as follows: the lowest G\ u 
state corresponds to spin \ and the lowest two H u states 
correspond to distinct spin | states. Above these is a 
group of three states with roughly the same lattice en- 
ergy: a t E w 0.33 (within errors). Again there are two 
possible interpretations. A.) The group consists of a spin 
\ state in G\ u that is accidentally degenerate with a spin 
| state, the latter having degenerate partner states in H u 
and G2u- B.) The group consists a spin | state having 
degenerate partner states in G\ u , H u and Cm- 

The pattern of low-lying physical states starts with 
A/"(§",1520), A(±~,1535) and 7V(§~,1700). These 
should correspond to distinct H u , G\ u and H u states on 
the lattice, in agreement with the three lowest negative- 
parity states in Fig. 4. The next two physical states are 

A/"(i",1680) and JV(f 1675), which essentially are de- 
generate. They should show up as degenerate G\ u , H u 
and Giu states on the lattice. This pattern of spins is 
consistent with interpretation A.) of the lattice states at 
lattice energy a t E s=a 0.33. The pattern of energies of the 
physical states has iV(| , 1675) a little lower in energy 
than ./V(| , 1700), but the lattice results at lattice spac- 
ing 0.1 fm place the spin | state above the spin | state. 
Study of the continuum limit of the lattice spectrum is 
required in order to resolve these issues. 

Figure 5 shows the pattern of lowest-lying energies for 
each irrep with the A states included. Energies have been 
converted to physical units but the purpose of this graph 
is to display the pattern of level orderings for low lying 
spin-parity channels with spins up to §, not to make a 
comparison of absolute energy values with experiment. 

The spin assignments in the negative parity channels 
obtained in our simulations correspond reasonably well 
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FIG. 5: The lowest energies obtained for each symmetry chan- 
nel of I = | and / = | baryons are shown in columns 2 and 
4 and the experimental levels are shown in columns 1 and 3. 
The 24 3 x 64 lattice data are shown and the left two columns 
show the energies of positive-parity states and the right two 
columns show the energies of negative parity states. The ver- 
tical axis shows energy in MeV. 



sponds, is slightly above, and finally the A(| ) is high- 
est. One result that does not correspond with nature is 
the relative order of the A(| ) and the (A| ) states, al- 
though the error bars overlap and the energy splitting of 
the mean values is only 2%. Considering these six states, 
the mean energies of iV(| ) and iV(| ) are less than 
the other four. This agrees with what is predicted by 
the quark model with the contact spin-spin interaction 
between quarks treated as a perturbation. 

In the positive-parity channels, with the exception of 
the Roper resonance, our spin assignments also follow 
the same orderings as occur in nature: the N(^ + ) is the 
lowest, the A(| + ) is next, the A(| + ) is next, for which 
H g and G2 g are degenerate, the A(i + ) is next up, and 

finally the A(| + ) is highest. The analysis of Isgur and 
Karl based on the quark model yields a nearly degenerate 
set of states corresponding to A(± + , 1710), iV(| + , 1720), 

and A(| + ,1680). Their calculation explains well the 
small energy splitting, ~ 40 MeV, between these states. 
Our results yield a fair agreement with this result. En- 
ergies of the first excited state of I = ^,Gi g , the first 
excited state of / = i, H g , and the lowest / = i, G2 g are 
0.398(6), 0.395(7), and 0.362(17), respectively. Closeness 
of first two energies above is remarkable, but the third en- 
ergy is significantly different. Note that our assignments 
cause the lowest A^(| + ) to correspond to the first excited 



with the physical spectrum, i.e., the A(^ ) is the low- state q{ j = l Hg because the lowest state corresponds 
est, the N(t£ ) is slightly above, the A(i ) is next, the |- J p 
A(| ) is next, the 7V(| ), to which / = A,G2 U corre- 



. Although our energy splittings at this 
quark mass can be much larger than those in the physical 
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spectrum, if our spin assignments are correct the order- 
ing of energies of the lightest states in a given channel 
reproduces the physical ordering. 

On the other hand, comparison of the order of en- 
ergy levels with different parities provides different con- 
clusions. For instance, the A(| ) energy should appear 

between AT(§ ) and AT(§ ) energies (here only the lowest 
energies are considered), but it docs not. Further study 
of the continuum limit is required in order to confirm the 
spin assignments. 

VII. SUMMARY AND OUTLOOK 

This paper presents spectra for isospin 1 = \ and 
7=| states based on lattice QCD simulations using 
the quenched approximation, anisotropic lattices with 
a s /a t — 3 and a pion mass of 490 MeV. Sets of three- 
quark interpolating fields are used to construct matri- 
ces of correlation functions and the variational method 
is used to extract energies. Smeared quark and gluon 
fields are used to diminish the coupling to short wave- 
length fluctuations of lattice QCD. The result is that 
diagonalizations are able to extract good signals for low- 
lying states. We have obtained as many as 17 energies 
for I = i states and 10 energies for / = | baryon states, 
including various spin-parity channels. The variational 
method determines the best linear combinations of basis 
operators for energy eigenstates. For each obtained en- 
ergy eigenstate, we paid close attention to the stability 
of the eigenvector components with respect to time and 
with respect to the two different volumes. 

Because scattering states are expected to be present 
in the spectrum, we have calculated spectral weights for 
all states on two volumes. However, the weights do not 
provide evidence for two-particle states in our spectra. 
The first positive-parity excited state in the I = \,G\ g 



channel is found to have energy significantly higher than 
the experimentally-known mass of the Roper resonance. 
Better ways are being developed to identify scattering 
states and resonances, such as the use of operators that 
should have large couplings to the tt — N states and the 
use of several pion masses in order to better detect s-wave 
wN states. 

Because the minimum spin that is contained in the 
irrep is |, we have found strong evidence for spin | or 
higher for both parities in our spectra for I = i. We 
also have found evidence for degenerate partner states 
corresponding to the subduction of spin | or higher to 
the octahedral irreps. 

The lattice results for the I = \ low-lying excited 
states provide the correct number and pattern of octahe- 
dral states for the subduction of the spins of the low-lying 
physical states, with the exception of the Roper state. 
When the A states are included, the overall pattern of 
lattice results for the lowest energy in each channel is 
similar to the pattern observed in nature. 

The results shown here and in Ref. [11] represent a first 
glimpse of the pattern of nucleon and A excitations as 
predicted by QCD, but they are based on the quenched 
approximation, a 490 MeV pion mass and an action that 
may have significant discretization effects. Similar cal- 
culations in full QCD are under way with an improved 
action, several pion masses, several lattice spacings and 
much larger sets of baryon interpolating operators. 
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